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Abstract 

We establish a uniform dimensional result for normally reflected Brownian motion 
(RBM) in a large class of non-smooth domains. Exact Hausdorff dimensions for the 
boundary occupation time and the boundary trace of RBM are given. Extensions to 
stable-like jump processes and to symmetric reflecting diffusions are also mentioned. 
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1 Introduction 



Let n > 2 and D C R n be a domain (connected open set) with compact closure. Consider 
a reflecting Brownian motion (RBM in abbreviation) X in D. Heuristically, RBM in D is 
a continuous Markov process X taking values in D that behaves like a Brownian motion 
in R" when X t G D and is instantaneously pushed back along the inward normal direction 
when Xt G dD. RBM on smooth domains can be constructed in various ways, including 
the deterministic Skorokhod problem method, stochastic differential equation with boundary 
condition, martingale problem methods, etc. see the Introduction of [Bj. When D is non- 
smooth, all the above mentioned methods cease to work. On non-smooth domains, the most 
effective way to construct RBM is to use the Dirichlet form method, which will be recalled 
in Section The RBM constructed through Dirichlet form coincides with all the other 
standard definitions in smooth domains. Using the Dirichlet form approach, it can be shown 
that, when D is a simply connected planar domain, RBM X in D is the time change of the 
conformal image of RBM in a unit disc. 

This paper is concerned with Hausdorff dimensions of various random sets associated with 
RBM. The study of Hausdorff dimensions of random sets associated with Brownian motion, 
stable processes and more generally Levy processes together with their fractal structures has 
been an active research area in the last 40 years. See Xiao [SHI for a recent survey on this 
subject. However, to the authors' knowledge, this is the first time that such a study has 
been conducted for RBM in Euclidean domains. 

The results obtained in this paper hold for a large class of non-smooth domains. In order 
to convey our results as transparently as possible, in this introductory section we confine 
ourselves to a special case of the more general results established in this paper by assuming 
that D is a bounded uniform domain. The following definition is taken from Vaisala [29J, 
where various equivalent definitions are discussed. 

Definition 1.1 A domain D C R n is called uniform if there exists a constant C such that 
for every x,y G D there is a rectifiable curve 7 joining x and y in D with length^) < C\x — y\ 
and moreover min {\x — z\, \z — y\)} < Cdist(z, dD) for all points z G 7. Here dist(z, dD) 
is the Euclidean distance between point z and the set dD. 

For example, the classical van Koch snowflake domain in the conformal mapping theory is 
a uniform domain in R 2 . The uniform domain is also called (e, 00) -domain in the terminology 
of Jones Note that every Lipschitz domain is uniform, and every non-tang entially 

accessible domain defined by Jerison and Kenig in JH] is a uniform domain (see (3.4) of [18J). 
However, the boundary of a uniform domain can be highly nonrectifiable and, in general, no 
regularity of its boundary can be inferred (besides the easy fact that the Hausdorff dimension 
of the boundary is strictly less than n). For any a G [n — 1, n), one can construct a uniform 
domain D C R n such that 7i a (U H dD) > for any open set U satisfying U D dD 7^ 0. Here 
7i a denotes the a-dimensional Hausdorff measure in R n . 
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In this paper, we first extend Kaufman's uniform doubling dimension result |20j for planar 
Brownian motion to RBMs in bounded domains. It follows as a special case of Theorem 13.11 
below that for RBM X in a bounded uniform domain D C R", 

P x (d\m. H X(E) = 2 dim H E for all Borel sets E C R+) = 1 for every x G D. (1.1) 

Here X(E)(u) := {X t (u) : t G E} denotes the range of E under RBM X and dim# denotes 
the Hausdorff dimension. Such a result is called a uniform dimensional result because the 
exceptional set in (jl.l|) is independent of the Borel time sets E C R + . This dependence is 
important when one wants to extract information on X(E) by observing E only while E 
itself is random, for example, when E is the boundary occupation time set of X. 

We next study the occupation time of X on the boundary dD, 

S(co) = {t>0: X t (oj) G 3D}. 
Corollary 14.31 below implies that 

d imH 5M = i- "- di ; n " ac (i.2) 

P^-almost surely for every x G D. Note that for any Euclidean domain D C R ra , dim^ dD G 
[n — 1, n]. This together with (jl.ljl implies that the Hausdorff dimension of the boundary 
trace for RBM in D is 

dim H (X[0, oo) n dD) = 2 + dim H dD - n 

P^-almost surely for every x G D. In particular, for planar uniform domains (such as the 
van Koch snowflake) it follows that P^-almost surely 

dim^(X[0, oo) n 3D) = dim H 3D 

for every x G D. This is in contrast to Makarov's celebrated result about the support of 
harmonic measure: There is a subset A of 3D with dim^ A — 1 such that the first intersection 
of X with 3D is almost surely contained in A. In Section El an example of a (non-uniform) 
domain D is given where the Hausdorff dimension of the boundary trace of RBM in D differs 
from dimiidD. Thus some assumption about the regularity of the domain is necessary for 
our results to hold. 

The remainder of this paper is organized as follows. In Section 2, we recall the defi- 
nitions of RBM and extension domains. In Section 3, we show that the sample paths of 
RBM have the same degree of Holder continuity as Brownian motion, by using the Lyons- 
Zheng's forward-backward martingale decomposition of RBM. We then establish the uniform 
dimensional results for RBM on a large class of non-smooth Euclidean domains. For this we 
derive the two-sided heat kernel estimates for RBM and use a capacitary argument. In the 
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case of dimension n — 2, the logarithmic capacity is not good for our approach so we first 
subordinate the RBM, then establish dimensional results for the subordinated stable-like 
process and lastly transfer these results back to RBM. In this procedure, we use the recently 
obtained two-sided heat kernel estimates in Chen and Kumagai [Zj for stable-like processes 
on D. In Section 4, we apply a recent result in Bogdan, Burdzy and Chen [3] to get the 
exact capacitary dimension of the boundary occupation time of RBM and therefore establish 
its Hausdorff dimension. All these results are combined in Section 5 to get the Hausdorff 
dimension for the boundary trace of RBM. Several examples are given to illustrate the main 
results of this paper as well as an example that shows that certain regularity assumptions 
on the domain D are needed for the main results of this paper to hold. Extensions of the 
main results in this paper to stable-like jump processes are mentioned in Remarks 13. 1(JI 14*741 
andO 

We point out here that the approach used in this paper is quite robust that applies not 
only to RBM but also to symmetric reflecting diffusions as well. These symmetric reflecting 
diffusions have infinitesimal generators of divergence form: 




where A(x) = (aij(x))i<i,j< n is a symmetric n x n matrix-valued measurable function that 
is uniformly elliptic and bounded. See Chen [B] for relevant informations on symmetric 
reflecting diffusions on bounded domains, and Stroock [28J for two-sided heat kernel estimates 
for symmetric diffusions on R™. 

2 Extension domains and RBM 

Throughout this paper, n > 2 is an integer. For a domain D C R n , let W 1,2 (D) := {/ G 
L 2 (D,dx) : V/ G L 2 (D,dx)}, equipped with the Sobolev norm ||/|| 1)2 := ||/|| a + ||V/|| 2 , 

1 /2 

where ||/||2 : = (f D f(x) 2 dx) . In this paper we will be concerned with reflecting Brownian 
motion on domains D C R™ whose boundary has zero Lebesgue measure and have the 
following W /1,2 -extension property: there is a bounded linear extension operator 

T : W 1,2 (D) -> W l ' 2 (R n ) such that Tf = f a.e. on D for / e W h2 (D). (2.1) 

The (e, 5)-domains (also called locally uniform domains) introduced by P. W. Jones [T7| 
have this property. A domain is (s,S), if there are s G (0, 00) and 5 G (0, 00] such that 
the conditions in Definition 11.11 for uniform domains hold for \x — y\ < 5 with C = 1/e. 
Clearly, an (e, oo)-domain is just a uniform domain. For later use, we record the following 
observation as a lemma. 

Lemma 2.1 If D C R n is an (e, 5) -domain, then so is D x [0, 1]. 
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Proof. If (x,t), (y, s) E Dx [0, 1] with \{x,t) — (y, s)\ < min(l, 5), then either \x — y\ > \t — s\ 
or \x — y\ < \t — s\. In the first case, we can easily lift the rectifiable curve 7 from x to y in D 
in the definition of (e, 5)-domain for D to a rectifiable curve from (x, t) to (y, s) that satisfies 
the (e, <5)-condition for D x [0, 1], with possibly different values of e and 5. In the second case, 
we construct a point z G D with \x — z\ = 2\t — s\ such that dist(z, dD)) > \t — s\/C, and 
then apply the first case to construct rectifiable curves joining (x,t) with (z, (t + s)/2) and 
(z, (t + s)/2) with (y, s). The piecing together curve satisfies the condition for (e, £)-domain 
for D x [0, 1], with possibly different values of e and 5. □ 

The importance of (e, #)-domains is further illustrated by the following. It is shown by 
Jones that a finitely connected domain is a H^ 1,2 -extension domain if and only if it is 
an (e, <5)-domain. More information on domains that have the W^ 1,2 -extension property (|2.1j) 
can be found in Herron and Koskela [To] . 

For f,ge W^iD), define 

£(f,g) :=^J D Vf-Vgdx, 

and £i(f,g) := £(f,g) + f D f(x)g(x) dx. We say the Dirichlet form {£ , W 1,2 (D)) is regular 
on D if W 1 ' 2 (D)nC(D) is dense both in (W^ 2 (D), £\ /2 ) and in (C(D), \\ ■ (U). It is known 
(cf. Theorem 2 on p. 14 in Maz'ja [24]) that (S, W 1,2 (D)) is a regular Dirichlet form on D if 
D has continuous boundary. Clearly (8, W l,2 (R n )) is a regular Dirichlet form on R ra . So if 
D is a iy 1,2 -extension domain, (S, W 1,2 (D)) is regular on D. 

When (£, W 1,2 (D)) is regular on D, there is a strong Markov process X associated 
with {£, W 1,2 {D)) ) having continuous sample paths on D; one can construct a consistent 
Markovian family of distributions for the process starting from every point in D except 
possibly for a subset iV of dD having zero capacity (see Chen [U]). When D is smooth, the 
exceptional set iV can be taken to be the empty set. We will show in Section 3 (see ()3.6|) 
and Lemma f3.7l below) that the exceptional set N in fact can be taken to be the empty set 
for a large class of nonsmooth domains including (e, 5)-domains. Thus constructed process 
X is the reflecting Brownian motion on D in the sense that this definition agrees with all 
other standard definitions in smooth domains. As it is remarked in Bass, Burdzy and Chen 
[2] (see the second paragraph in the proof of Theorem 5.7 there), such reflecting Brownian 
is conformally invariant on planar domains. 

Let S(u) := {t > : X t {u) G dD} be the occupation time of X on the boundary dD 
and R(u) := X[0, 00) (ou) R dD the trace of X on the boundary dD. 

Recall that for any increasing function h on [0, 1) with h(0) = 0, one can define a 
Hausdorff measure Hh with respect to the gauge h in the following way (see, e.g., p. 132 of 
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[U). For E C R n , 

{oo oo 
h(rk) '■ E C B(xk, rfc) for some G R n with sup < £ > . 
fc=l fc=l l<fc<oo J 

(2-2) 

When /i(r) = r a for some a > 0, the Hausdorff measure Tih will be denoted as H, a . 

Definition 2.2 A Borel set T C R n called an n-set if there exists a positive constant 
c > suc/i i/iai 

m(T fl r)) > cr n for all x G T and < r < 1, 

where m denotes the Lebesgue measure in R n . 

Note that if T C R n is an n-set, then by Proposition 1 in Chapter VIII of so is 
its Euclidean closure T and T \ T has zero Lebesgue measure in R n . It is known that any 
(e, <5)-domain in R n is an n-set (see Example 4 on page 30 of |19j). 

3 Uniform Dimensional Results for RBM 

In this section, we will extend Kaufman's uniform dimensional result for planar Brownian 
motion to reflecting Brownian motions in n-dimensional bounded domains. This is perhaps 
the first time such a uniform dimensional result has been established for a recurrent Markov 
process that does not have the property of independent increments. 

Theorem 3.1 Suppose that D C R" is a bounded domain whose boundary dD has zero 
Lebesgue measure such that either 



(1) it satisfies the extension property \2.1\) when n > 3; 



or 



(2) D C R 2 is a connected open 2- set such that the product domain D x [0, 1] satisfies the 
extension property \2.1\) . 

Then for every x G D, 

P x (dim H X(E) = 2dim H E for all Borel sets E C R+) = 1. (3.1) 

We will prove the upper bound first. 

Theorem 3.2 Let D C R n be a domain such that {S,W^ 2 {D)) is a regular Dirichlet form 
on D. Then 
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(i) for every x G D 



P x (dim H X(E) < 2 dim H E for all Borel sets E C R+) = 1; 



(3.2) 



(ii) if D C R n satisfies the condition of Theorem \8.1\ then h8. jjj) holds for every x G D. 

Proof. We first show that the sample paths of RBM have the same degree of Holder 
continuity as Brownian motion. Clearly, since RBM X t behaves like Brownian motion when 
X t G D, the sample paths of X can not be smoother than those of Brownian motion. 
By Lyons-Zheng's forward and backward martingale decomposition (see Theorem 5.7.1 of 
Fukushima, Oshima and Takeda [TT|). for any T > 0, 



where If is a martingale additive functional of X which is an n-dimensional Brownian 
motion, and tt is the time reversal operator of X at time T, i.e., X t (rT{uj)) = for 
each < t < T. Since X is symmetric under P m , P m is invariant under time reversal Tt- 
For a G (0, ~), define 



It follows from ()3.3j) that f a — 1 q.e. on D. Here q.e. is the abbreviation for quasi- 
everywhere; that is the above property holds holds for every point x in D except a set of 
zero capacity with respect to the RBM X. Since / is finely continuous, we have f a = l q.e. 
on D. This implies that for q.e. 

P x (dim H X(E) < or 1 dim H E for all Borel sets E C R+) = 1. 

Taking a sequence a n f 1/2, we have for q.e. x G D, 

P x (dim H X(E) < 2 dim H E for all Borel sets E C R+) = 1. 

Since X has density function p(t, x, y) for every x G D (see Fukushima PHI), the above holds 
for every x G D. If D C R n satisfies the condition of Theorem 13. 11 then by ()3.6|) and Lemma 
13. 71 below. X has transition density function p(t,x,y) for every x £ D. It follows that ()3.2j) 
holds for every x E D. □ 

Remark 3.3 In fact, Theorem I3.2IT J holds for RBM on any Euclidean domain, with the 
same proof. See Chen |H| for the construction of RBM and its forward-backward martingale 
decomposition on an arbitrary Euclidean domain. When D is a domain satisfies the condition 
of Theorem 13.11 ()3.2|) can also be established by using the heat kernel estimates obtained in 
(13. 6 j) and Lemma f3 . 71 together with a result initially due to Hawkes and Pruitt (see Lemma 
8.1 of Xiao (SO!. □ 




(W T - W T -t) o r T for all < t < T, P m -a.s., 



(3.3) 



fa(x) 



P x (t — ► Xt is a-H61der continuous on each finite interval) . 
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To prove the lower bound in Theorem 13. 1[ we treat the higher dimensional case (n > 3) 
and two-dimensional case separately. 

Let D be a domain that satisfies the condition of Theorem 13.11 Let N be a Borel set of 
D having zero capacity with respect to RBM X in D such that for every x G D\N, 

P x (there is some t > such that X t G N or X t _ G N) = 0. 

Such a set N is called a properly exceptional set of X and alway exists (see [TT]). When 
n > 3,by the classical Sobolev inequality on R" and (|2.1j) . 

||/|| 2 n/(n-2) <cy/£jjj) for / G ^(D). (3.4) 

Thus according to Varopoulos (see Theorem 2.4.2 in Davies (HJ) the reflecting Brownian 
motion X has density function p(t, x, y) such that 

e~ l p{t,x,y) < cr n/2 for t > and x,y G D \N. (3.5) 

Recall that it is assumed that 3D has zero Lebesgue measure. Using Davies' method together 
with an old idea of Nash advanced by Fabes and Stroock in j^j, it is now standard to deduce 
that (cf. Theorem 2.3 and Theorem 3.4 in [H| and Section 3 in jHj) p(t,x,y) is jointly Holder 
continuous on R + x D x D and for any k > there are positive constants ci, C2, C3 and C4 
such that 

ci r"/ 2 exp (-^^) < P(t, x, y) < c 3 r n ' 2 exp (Jl^Pj (3.6) 

for < t < k and x,y G D. So p(t, x, y) can be extended continuously to [0, 00) xDx D. This 
in particular implies that reflecting Brownian motion on a W /1 ' 2 -extension domain D C R n 
with n > 3 can be constructed as a strong Markov process that starts from every point in 
D(cf. HU). 

The following covering lemma was first proved for Levy processes in R n by Pruitt and 
Taylor [27|. It is extended to general Markov processes in Liu and Xiao [22J. 

Lemma 3.4 Let A(a) be a collection of cubes of side a G (0, 1] in R n with the property that 
the number of these cubes which intersect an arbitrary sphere of radius a in R™ is bounded 
by a constant K that is independent of a and of the sphere (this happens when the cubes are 
those in aZ n or when the cubes do not overlap too much). Let D C R n satisfy the condition 
of Theorem \3.1\ and let M(a,t) be the number of those cubes that are hit by RBM X in D 
before time t. Then there is a constant c = c(K,t) that depends only on K and t such that 

W [M(a,t)\ < c 

for every x G D. 



inf_ 

y&D 



p(s, y, z) dz ds 



JB(x,a/3) 



(3.7) 
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It follows from ()3.6|) that 

E x [M{a,t)}<ca- 2 . (3i 

Lemma 3.5 If U is a ball of radius a in R n and k > 0, then there is a constant c ■ 
c(D, k) > such that 

P x (X s G U for some s G [t, k}) < c (^f^ 
for every a < 1 and x G D. 

Proof. Let be the 1-equilibrium measure of X for U n D; that is, 

E*[e-^]= f G 1 (x,y) f x(dy), 
Jd 

where Ou := inf{t > : X t £ U} and Gi(x,y) := / °° e~ t p(t,x,y)dt. Note that 

P x (X s G C/ for some s G [t, fc]) < E x [P Xt (X s G U for some s G [0, k])] 

< e k E x [E Xt [e~ au ]] 



< e k J_{j e (s t) p(s,x,z)ds) fi(dx) 



< ce 2k J s- n/2 dsfi(U) 

= c r("- 2 )/ 2 Ca Pl (f/). 

It follows from the extension property (j2.1|) that Cap 1 (?7) < cCapf' (U), where Cap 1 (f/) 
and Capf^ (U) denote the 1-capacity of U for RBM X on D and Brownian motion in R n , 
respectively. Hence 

Ca Pl (£7) < ca n ~ 2 , 

which proves the lemma. □ . 

Theorem 3.6 Suppose that n > 3, that D C R n is a bounded domain satisfying the exten- 
sion property \2.1\) and that 3D has Lebesgue measure. Then for every x G D, 

P x (dim H X(E) > 2 dim H E for all Borel sets E C R+) = 1. (3.9) 

Proof. With Lemmas 13 .41 and 13 . 51 in hand, by the same argument as those for Lemma 3 and 
Theorem 1 in Hawkes [T3j . we have 

P x (dim H I(£n[0,fc])>2 dim H (E n [0, k}) for all Borel sets E C R + ) = 1 

for every k > and x G D. Letting k — > oo, we see that fl3.9j) holds for every x G D. □ 

So far we have proved Theorem 13. lf l ) . In the remainder of this section, we will deal with 
the two-dimensional case. 
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Lemma 3.7 Suppose D C R 2 is a bounded domain whose boundary 3D has zero Lebesgue 
measure such that D x [0, 1] is a W 1,2 -extension domain in R 3 . Then reflecting Brownian 
motion on D exists as a strong Markov process on D starting from every point in D. Fur- 
thermore its transition density function p(t, x, y) is jointly Holder continuous on R + x Dx D 
and has Gaussian estimates \3. 6}) with n = 2 there. 

Proof. Under the assumptions of the lemma, D is a W 1,2 -extension domain in R 2 so RBM 
X on D exists as a strong Markov process on D starting from every point in D except on a 
properly exceptional set N of X. Note that RBM Y in D x [0, 1] can be obtained from X 
by running an independent RBM in the unit interval along the ^-direction. By ()3.5|) . RBM 
Y on D x [0, 1] has transition density function p(t, x, y) and 

e~ l p(t, x, y) < ct~ 3/2 for every t > and x, y e (D\N) x [0, 1]. 

It follows that X has density function p(t, x, y) and that 

e~*p(i, x, y) < ct^ 1 for every t > and x, y G D \ N. 

Joint Holder continuity of p(t, x, y) and its two-sided estimate (|3.6|) (with n = 2) now follows 
from a similar argument as that in Section 3 in Fabes and Stroock [9J and that for Theorems 
2.3 and 3.4 in Bass and Hsu 0. So the RBM X can be refined to start from every point in 
D. □ 

Remark 3.8 By Lemma 12.11 and the remark at the end of Section 2, any bounded planar 
(e, 5)-domain satisfies the condition of Lemma 13.71 Combining this with a result of Jones 
mentioned at the paragraph following Lemma I2~T| we see that any bounded finitely connected 
H /1 ' 2 -extension domain in R 2 satisfies the condition of Lemma f3. 71 □ 

Theorem 3.9 Suppose that D C R 2 is a connected bounded open 2-set that satisfies the 
condition of Lemma \3. 7| Then for every x G D, \3. .9)) holds. 

Proof. Let X be a reflecting Brownian motion in D with sample space (Q, P x ,x G D). By 
Lemma l3~7l the transition density function has the estimate ()3.6|) with n = 2. For s G (0, 1), 
let £ be a R + -valued process on a probability space (Q', P') with independent increments 
such that 

E [e~ A€t ] = e~ txa for any A > 

We assume that X and £ are independent and live on a common probability space. This can 
be achieved on the product space Q x Q with product measures {P x (g> P', x G D}. Define 
Z t = X^ t ; that is, Z is the s-subordinate of X. 
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It follows from Kumagai [21] that Z is a special case of the stable-like processes on D 
studied in Chen and Kumagai [7j. Hence by Theorem 1.1 of [7j the transition density function 
h(t, x, y) of Y has the property that 

e~*/i(t, x, y) < ct~ 1/s for t > and x,y E D 

and 

h(t,x,y) « min <^ — — I for ijgI) (3.10) 

[ |x — J 

on any finite time intervals. Here for two function f,g, f w g means there are two positive 
constants c\ < C2 such that c±g < / < C2#. Furthermore, it follows from Theorem 2.5(1) and 
its proof in Bogdan, Burdzy and Chen [3j that the capacity induced by Z on D is controlled 
by the capacity induced by the symmetric (2s)-stable process in R 2 (cf. also Lemma 3.1 of 
Fukushima and Uemura J2])- Thus by a similar argument as that for Lemma f3.5l above, we 
have 



P x ® P' (Z r G U for some r G [t, k]) < c 




2-2s 



where U is a ball of radius a in R n , k > 0, and c = c(D, k) > is a constant that depends 
only on s and k. Using this and applying Lemma I3~^l to the process Z, the same reasoning 
as that in the proof of Theorem 13.61 implies that for every x G D, 

P x ® P' (dim^ Z(E) > 2s dim H E for all Borel sets E C R+) = 1. (3.11) 

For a Borel set E C R+, define C E = {t > : & G £}. Clearly dim H X(E) > dim H Z(C B ), 
and the random number dim^C^ depends only on the subordinator £ and E, which is 
independent of X and therefore of dini^ X(E). Thus by ()3.11j) and Fubini's theorem, for 
every x G D, there is Q C Q with P x (n ) = 1 such that for every u & Q and for every 
Borel measurable set E C R + , 

dim H X(E)(u)> 2s dim H C E (u') for P'-a.s. u/ G f2'. (3.12) 

Thus for every u; G fi and every Borel measurable set C R + , 

dim H X(E)(u) > 2s sup{6 : P^dim^ C B > b) > 0}, 

which by (3.8) of Pruitt j2H] is no less than 2(s + dimjy — 1). Hence we have shown that 
for every x G D, 

P x {dim H X(E) > 2(s + dim H E - 1) for all Borel sets E C R+) = 1. 
Letting s f 1 proves the theorem. □ 

Theorems 13.21 13.61 amd 13.91 imply Theorem 13.11 
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Remark 3.10 Let T C R n be a closed set such that there are constants d G (0, n] and 
c 2 > C\ > so that for every x G T, 



H d (B(x, r) n T) > Cir d for < r < 1 and H d (B(x, r) n T) < c 2 r d for r > 



(such a set is called a d-set). Let Y be an a-stable-like process on V (see Chen and Kumagai 
[7j for the definition), where < a < 2 and a < d. It is shown in [7] that F has a jointly 
Holder continuous transition density function p(t, x, y) on [0, oo) xfxT and that 



on [0, k] x T x T for every k > 0. Using a similar argument as above, one can show that 

P x (dim# Y(E) > a dim^ E for all Borel sets E C R + ) = 1 for every x ET. 

As it is observed by Xiao in [SU], applying a result initially due to Hawkes and Pruitt (see 
Lemma 8.1 of [SO]), one gets 

P x (dinifl- Y (E) < a dim^- E for all Borel sets E C R + ) = 1 for every x G T. 

Hence we have for < a < 2 and a < d, 

P x (dim H Y(E) = a dim H E for all Borel sets E C R+) = 1 for every x G T. (3.13) 

Taking E = R + , one can see clearly that the above can not be true when a > d. □ 

4 Boundary Occupation Time 

Theorem 4.1 Suppose that D is a bounded open n-set satisfying condition \2.1\) . Then for 



(i) dim# S(uj) < 7 P x -a.s. for every x G D if 3D is polar for the symmetric 2(1 — ^-stable 

process in R n . 

(ii) dini^ S(uj) > 7 P x -a.s. for every x G D if 3D is not polar for the symmetric 2(1 — 7)- 
stable process in R n . 

Proof. For < s < 1, let Z be the s-subordinator that is independent of RBM X in 
D. Then Y = Xz t is a recurrent symmetric Markov process with density functions and its 
associated Dirichlet form has the property that its ^-norm is comparable to the Besov norm 
in Bl' 2 (D) (see Kumagai for a calculation). Hence by Theorem 2.5 in Bogdan, Burdzy 
and Chen j3] , dD is polar for Y if and only if it is polar for the symmetric 2s-stable process 




7 G(0,1) 
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in R n . Note that according to Frostman the Hausdorff dimension is the same as the Riesz 
capacitary dimension, while Orey showed that one- dimensional a-stable processes share the 
same polar sets (see Lemma 2 in Hawkes [Oj). Therefore we have for any Borel set E C R+, 

dim// E = sup {1 — s > : E is not polar for the s-subordinator Z} . 

So for every x G D, P x -a..s., 

dim// S(lu) 



sup {1 — s > 
sup {1 — s > 
sup {1 — s > 



dD is not polar for the s-subordinator Z} 
dD is not polar for Y} 

dD is not polar for the symmetric (2s)-stable process in R n } . 



This proves the theorem. □ 

The above theorem implies that dim// S (uj) = 7 almost surely if and only if dD is polar 
for symmetric a-stable process in R n with a < 2(1 — 7) an d non-polar for symmetric a-stable 
process with a > 2(1 — 7). 

Remark 4.2 There is an intimate relationship between the Hausdorff measure Hh (see ()2.2|) 
for its definition) and the Riesz capacity Cap ra _ Q of order n — a, see Theorems 2.2.7, 5.1.9 
and 5.1.13 in PQ. Namely, TL n ~ a (A) < 00 implies that Cap n _ a (A) = 0. On the other hand 
if Cap ri _ Q ,(^4) = then 7i h (A) = for every h such that 

h is increasing on [0, 00) with h(0) = and / _ q dr < 00. (4-1) 

Jo r ™ 

In particular, C&p n _ a (A) = implies 7i x (A) = for any A > n — a. Later we will use the 
fact that for n = 2 and h(t) = 1/ log(l/t), 7ih(A) = implies cap A = 0, where cap A stands 
for the logarithmic capacity see Theorem 5.1.9 in . □ 



Recall that for any Euclidean domain D C R n , dim// dD e [n— 1, n\. Combining Remark 
14. 21 with Theorem 14.11 we have the following corollary. 

Corollary 4.3 Suppose that D C R n is a bounded open n-set satisfying condition \2.1\) . 
Then P x -a.s. dim// S (uj) — 1 — | (n — dim// dD) for every x e D. 

Proof. If for some d G [n — l,n), 7i d (dD n K m ) < 00 for an increasing sequence of Borel 
sets K m such that U™ =1 K m D dD, then by Remark l4~2} C&p d (dD) = and so dim h S(lj) < 
1 - ^ almost surely. If H d {dD) > for some d G [n — 1, n), then Cap d (<9-D) > and hence 
dim// S(to) > 1 — -^rp almost surely. Since 

dim// dD = inf {a > : H a (dD) = 0} = sup{a > : ft Q (0£>) = 00}, 

the conclusion of the corollary now follows. □ 
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Remark 4.4 Let Y be an a-stable-like process on an open n-set D C R n in the sense 
of Chen and Kumagai [?]. It includes as a special case the reflected a-stable process on 
D introduced in Bogdan, Burdzy and Chen jl]. A similar argument gives the Hausdorff 
dimension of the boundary occupation time for Y, which asserts that P^-almost surely 

dimfl- S Y (lu) = sup {1 — s > : dD is not polar for symmetric (as)-stable process in R n } 

for every x G D. Hence we have for every x G D, P^-a.s. 

dim H ^M=ma X [l- n - di ^^, o} . (4.2) 

Note that when dD has locally finite d-dimensional Hausdorff measure for d < n — a, it 
follows from that dD is a polar set for Y and therefore S Y (u) is the empty set almost 
surely. □ 



5 Boundary Trace of RBM 

Combining Theorem 13. II with Corollary 14.31 establishes the following result. 

Theorem 5.1 Let D C R n with n > 2 be a bounded connected n-set satisfying the condition 
of Theorem \S.l\ Then the Hausdorff dimension for the boundary trace of RBM X in D is 

dim H (X[0, oo)ndD) = 2 + dim H dD - n P x -a.s. 

for every x G D. 



Remark 5.2 Let D C R n be an open n-set and Y be an a-stable-like process on D in the 
sense of with a < n. Then it follows from Remarks 13 . 1 01 and 14 . 41 1 hat the boundary trace 
R Y {u) := Y[0, oo) fl dD of Y has 

dim// R Y (ui) = max {a + dim^ dD — n, 0} P^-a.s. 

for every x G D. □ 



Example 5.3 Let D C R n with n > 2 be a bounded Lipschitz domain. Then it is an n-set, 
satisfies the condition of Theorem 13.11 and dim# dD = n — 1. So the boundary occupation 
time set and boundary trace of RBM in D have Hausdorff dimensions 1/2 and 1, P^-a.s., 
respectively, for every x G D. 
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Example 5.4 Let D C R 2 be a van Koch snowflake domain. By Remark 13.81 it satisfies 
the condition of Theorem 13. If 2) and it is well known that 3D has Hausdorff dimension |^||. 
Hence by Theorem 15.11 the boundary occupation time set and the boundary trace of RBM 
in D have Hausdorff dimensions |j^|§ and P^-a.s., respectively, for every x G D. 

Example 5.5 Let U = D x (0,1) C R 3 , where D is the van Koch snowflake domain in 
R 2 . Then U satisfies the condition of Theorem I3.1f 1) and dim// 3U = 1 + Hence by 
Theorem 15.11 the boundary occupation time set and the boundary trace of RBM in U have 
Hausdorff dimensions and P x -a.s., respectively, for every x G U. Note that the 

reflecting Brownian motion in U is (X, Y), where X is the RBM in D and Y is the RBM in 
(0, 1). The above boundary trace result may be a bit surprising if one compares it with the 
range and the graph of a 1-dimensional Brownian motion B. It is known that the Hausdorff 
dimension for the range {B t (u) : t G R+} of B is 1 a.s., while the Hausdorff dimension for 
the graph {(t, B t {u)) : t G R+} of B is 3/2 a.s. 

Example 5.6 Generalizing Example 5.4, let D C R 2 be a simply connected domain whose 
boundary is a regular fractal in the sense of (the typical examples are from complex 
dynamics, namely components of the Fatou set of a hyperbolic rational function, and self- 
similar curves such as the van Koch snowflake). It is well-known that 3D has positive and 
finite (i-dimensional Hausdorff measure, where d = dimn 3D. By Theorem 15 . 1 1 1 he boundary 
trace of RBM in D has Hausdorff dimension d almost surely. We will now outline an 
alternative proof of this result, based on complex analytic methods. Let / : D — > D be 
a conformal map of the unit disc D onto D. Using the methods and results from [23J, 
Section 3 (in particular Proposition 3.2), one can show that there is a set A C <9D with 
< dim// A < 1 such that 

dim H f(B) = d for all B C A with dim H B = dim H A. (5.1) 

In fact, one can show that there is a G (0, 1) such that the set 

{ r-i|log(l-r)| J 

satisfies dim# f(A) = d, and that / is uniformly expanding on A, 

C-i\ x - yl 1 -^ < \f(x) - f(y)\ < C e \x - 

for every e > and all x, y G A. This easily implies that A satisfies (15. 1|) . By conformal 
invariance /~ X (X) is a time change of RBM in D. Therefore f~ x (X[0, oo) PI 3D) has the 
same law as the boundary trace of RBM in D. Because RBM in D can be constructed 
from 2-dimensional Brownian motion by reflecting the part outside D, and because Brow- 
nian motion intersects A in a set of dimension dim// A almost surely by [25J, we obtain 
dim// f- 1 (X[0, oo) n 3D) = dim H A a.s. Now (JSHJ) implies dim// (X[0, oo) n 3D) = d al- 
most surely. 
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Figure 1: A finite approximation to B (black squares) and D (white squares) 



Example 5.7 Our final example shows that some assumption on the regularity of D has 
to be made in order for Theorem 15. II to hold. We will construct a simply connected planar 
domain D such that 

dD = A U B, 

where A has a-finite length, B has positive area, and f~ l (B) has zero logarithmic capacity, 
where / : D — > D is a conformal map. Thus planar Brownian motion does not visit f~ l {B). 
Using conformal invariance as in Example 5.6, it follows that X[0, oo) C D U A and hence 
dim^ (X[0, oo) n dD) = 1 a.s., while dim H dD = 2. 

The idea is to construct D as a countable union of "smooth" domains, for instance 
squares, joined by narrow "corridors". The limit set of the squares is B, and the width of 
the corridors can be arranged so that cap f~ x (B) = 0, where cap(A) denotes logarithmic 
capacity of a set A C R 2 . To give a rigorous description, it is easier to begin with the set B: 
Let B C R 2 be a totally disconnected compact set of positive area in the unit ball £>(0, 1). 
Consider a decomposition 

B?\B = \JS n 

n 

by closed squares S n with pairwise disjoint interiors, for instance the Whitney decomposition 
of R 2 \ B. It is easy to find line segments Ik such that each Ik is contained in some edge of 
some S n and such that 

D:= l{JS n U \JlA n£?(0, 2) 

\ n k / 

is connected and simply connected, see Figure 1. Then dD = A U B where A C (U n dS n ) U 
dB(0, 2) has u-finite length. Let Ik denote the length of Ik and / : D — ► D a conformal 
map. Since dD is locally connected, / extends continuously to D. Given any sequence 
tk > 0, it is easy to choose the Ik such that diam/ _1 (/ fc ) < tk for all k (for instance using 
Beurling's projection theorem, which gives diam/ _1 (/fc) < C\/lk). Denote Jk the arc on 
dD that is separated from by /~ 1 (/fc), so that diamJ fc < 2efc. If x G /^(B), then f[0,x) 
passes through infinitely many Ik. Therefore f~ l (B) C U n >kJ n for all k. By choosing 
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such that E^iVMVen) < oo we get H h (f~ x (B)) = with h(t) = l/log(l/t). Now 
cap f~ l (B) = follows from Remark 14.21 

Acknowledgements. We thank Y. Xiao for very helpful comments on a preliminary version 
of this paper. 
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